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Charged particle decay of hot and rotating ®®Mo nuclei in fusion-evaporation reactions 
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A study of fusion-evaporation and (partly) fusion-fission channels for the ®®Mo compound nucleus, 
produced at different excitation energies in the reaction ^^Ti A ®®Ca at 300, 450 and 600 MeV beam 
energies, is presented. Fusion-evaporation and fusion-fission cross sections have been extracted and 
compared with the existing systematics. Experimental data concerning light charged particles have 
been compared with the prediction of the statistical model in its implementation in the Gemini-|--|- 
code, well suited even for high spin systems, in order to tune the main model parameters in a 
mass region not abundantly covered by exclusive experimental data. Multiplicities for light charged 
particles emitted in fusion evaporation events are also presented. Some discrepancies with respect 
to the prediction of the statistical model have been found for forward emitted a-particles; they may 
be due both to pre-equilibrium emission and to reaction channels (such as Deep Inelastic Collisions, 
QuasiFission/QuasiFusion) different from the compound nucleus formation. 

PACS numbers: 25.70.Gh,25.70.Jj,24.60.Dr 


I. INTRODUCTION 

The study of the de-excitation of compound nuclei 
dates back to the beginning of the modern nuclear 
physics, but it has found renewed interest in recent years 
(for example, [T] for medium mass systems and [5]-[^ 
for very light systems) thanks to the improvement of 
the experimental techniques, allowing for a cleaner se¬ 
lection of exclusive decay channels. In this way it is pos¬ 
sible to put stringent constraints on the ingredients of 
the statistical model, commonly used to describe the de¬ 
cay of a compound nucleus at excitation energies below 
3 MeV/nucleon. The hypothesis behind the use of this 
model is that compound nuclei are equilibrated systems, 
whose decay is independent of their previous history (i.e. 
the way in which they were formed), but depends on the 
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statistical competition among all the open channels. The 
decay width of each channel depends on the phase space 
available in the final state and on the amplitude of the 
transition matrix. The tuning of the parameters of the 
statistical model on the basis of the experimental data 
obtained from compound nucleus reactions is important 
not only for the topic of low energy fusion reactions, but 
also for other classes of nuclear processes. In fact many 
kinds of reactions, from low (5-10 MeV/nucleon) to high 
(>100 MeV/nucleon) beam energies, such as deep inelas¬ 
tic collisions, multifragmentation, spallation, etc., can be 
theoretically described by means of models of the inter¬ 
action phase (e.g. dynamical models such as BNV [0], 
AMD [7], etc.) followed by a statistical afterburner for 
the hot products. Since in many cases the effects related 
to the nuclear interaction phase under investigation are 
weak and blurred by the afterburner, it is important to 
describe the statistical decay as well as possible in or¬ 
der to reduce the uncertainties on the dynamical stage. 
A good control of the decay process is also mandatory 
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in GDR (Giant Dipole Resonance) or DDR (Dynamical 
Dipole Resonance) studies because the collective excita¬ 
tion features are extracted from continuous spectra dom¬ 
inated by the statistical gamma rays from the source, 
which must be therefore correctly described by models. 

The present work deals with the investigation of fusion- 
evaporation and fusion-fission channels for the ®®Mo com¬ 
pound nucleus at three different excitation energies, pro¬ 
duced in the reaction '*®Ti-|-"‘°Ga at 300, 450 and 600 MeV 
beam energy. In particular, we compared the experimen¬ 
tal energy spectra and multiplicities of light charged par¬ 
ticles with the prediction of Gemini-|— I- code [8], a widely 
used implementation of the statistical model, in order to 
constrain the relevant parameters in a mass-energy-spin 
region in which not many experimental data of exclusive 
type, with a clean selection of the fusion channel, are 
available. Moreover, in this work we have also extracted 
the experimental cross sections for the fusion-evaporation 
and the fusion-fission channels, obtaining values compat¬ 
ible with the existing systematics [S]. 

This experiment is part of a measurement campaign 
aimed at the investigation of hot, rotating ®®Mo nucleus. 
The results concerning GDR gamma decay are presented 
in mm- The ®®Mo nucleus, as other medium mass nu¬ 
clei, is a good candidate for such kind of study because 
it presents a significant fission barrier even at high rota¬ 
tional energy, and it consequently has low fission proba¬ 
bility also at high spin values [12]. The use of GeminiH— h 
as statistical code to be compared with our experimental 
data is thus a good choice, since this code is particularly 
well suited for the decay of high spin nuclei. 

Among the data available in literature for the mass re¬ 
gion A < 100 we can cite, for example, [T3|, where the 
inclusive decay of ®’^Ga at 91 MeV excitation energy has 
been investigated. Proton and a-particle spectra at some 
polar angles have been compared with the prediction of 
a statistical code; a good agreement for a-particle cen¬ 
ter of mass energy spectra is obtained if a deformation 
of the source is taken into account. On the contrary, 
proton energy spectra are not well reproduced. Instead, 
the angular distributions fit quite well with the predic¬ 
tion of the code for both protons and a-particles, when 
a deformed source is used. 

An inclusive study of two different reactions producing 
the same compound nucleus ®^Ga at 1.9 MeV/nucleon ex¬ 
citation energy is presented in [14] . The authors found 
that it is not possible to describe both the energy spectra 
and the angular distributions of a-particles by means of 
a statistical code with a unique set of parameters; more¬ 
over, proton energy spectra require reduced evaporation 
barriers to be well reproduced. 

In [15] proton and a energy spectra in coincidence with 
evaporation residues have been measured for the reac¬ 
tion -I- ’^"‘Ge at different beam energies in the range 
5-13 MeV/nucleon. From the comparison with statisti¬ 
cal model predictions, the reverse level density parameter 
K = ^ is found to be substantially independent of the 
excitation energy of the system in the mass region around 


A = 100. The dependence of K on the angular momen¬ 
tum has been investigated in [T] for the compound nuclei 
®^Tc and ®^Zn at 36 MeV excitation energy, comparing 
the energy spectra of neutrons, protons and a-particles 
with the prediction of the Cacarizo code (Monte Garlo 
version of the Cascade statistical model code (l^); the 
angular momentum selection was performed by means of 
the 7 multiplicity. The authors found that the deforma¬ 
tion of the source (taken into account by means of the 
deformability parameters included in the code) has no 
effect on neutrons and protons energy spectra, while it 
significantly modihes those of a-particles for the lighter 
compound nucleus. The K values extracted from pro¬ 
ton, neutron and a-particle spectra are different, but they 
all show a decrease when the angular momentum of the 
source increases (in the range 13-22 h). 

In [17] a energy spectra detected in coincidence with 
Evaporation Residues (ER) from ®®Ni CN produced in 
the symmetric reaction ^®Si -I- ^®Si at 112 MeV beam en¬ 
ergy are compared with the prediction of Cacarizo. 
The experimental energy spectra have been well repro¬ 
duced by introducing a high degree of deformation (cor¬ 
responding to a quadrupole deformation parameter equal 
to 0.5) which modifies the effective moment of inertia of 
nuclei, thus lowering the yrast line. In this way the decay 
chain becomes longer and a-particles are emitted later in 
the cascade, with smaller average kinetic energy. The au¬ 
thors claimed also that deformation at high spin values 
is particularly favoured if the entrance channel is mass 
symmetric. Note that the spin range explored by this 
system is lower than in the present ®®Mo case; in fact the 
critical angular momentum for fusion was 34 h in m, 
while it is 64 h in our case. 

In [TH] a detailed study of the transmission coefficients 
of the statistical model, based on the energy spectra of 
hydrogen isotopes and a-particles evaporated by com¬ 
pound nuclei of ®®Ru at 1.2 MeV/nucleon excitation en¬ 
ergy, is presented, a-particle spectra are found to be 
in good agreement with the prediction of the statistical 
model (Cascade code) if reduced barriers with respect 
to the prescription of the optical model are used. More¬ 
over, the level density has to be enhanced, indicating the 
onset of deformations at high spins (kept into account 
by lowering the yrast line as done in m)- For protons, 
on the contrary, the authors explain the unsatisfactory 
agreement between experimental and simulated energy 
spectra with the need of taking into account dynamical 
effects. 

In [in] and [20] spectra of a-particles emitted from ^®Ti 
GN with 85 MeV excitation energy and critical angular 
momentum of 35 h are discussed and compared with the 
prediction of the already cited Cacarizo code. Various 
parametrizations of the yrast line implying different de¬ 
grees of deformation of the CN have been tested, obtain¬ 
ing reasonable agreement with the experimental data. In 
particular, energy spectra of a-particles detected in co¬ 
incidence with a residue oi Z = 20 (corresponding to 
chains in which only one a-particle is emitted) can be 
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reproduced only if an extremely high value of the de¬ 
formation parameter is used. The authors interpret this 
fact as due either to the presence of a dynamical hyper- 
deformed state or to the pre-equilibrium emission from 
a dinuclear system not completely fused. This study is 
coupled to the investigation of the GDR, confirming the 
presence of a strongly deformed structure. 

In summary, all these studies demonstrate that the de¬ 
cay features of light systems with high spin (up to ^ 60 Ti) 
are a difficult task for available statistical codes. In fact, 
a single parametrization is often not able to reproduce 
both proton and a-particle spectra. Moreover, in many 
cases a high degree of deformation is necessary to obtain 
reasonable results. We anticipate that also our results 
show that the statistical model with different assump¬ 
tions of the source parameters is not able to reproduce 
all the characteristics of the emitted particles, in particu¬ 
lar of the a particles, with increasing bombarding energy. 


II. EXPERIMENTAL SETUP AND ION 
IDENTIFICATION 

®®Mo compound nuclei have been produced at 
three different excitation energies (1.4, 2.2 and 

3.0 MeV/nucleon) by means of fusion reactions per¬ 
formed at the INFN Laboratori Nazionali di Legnaro 
(LNL, Italy). Pulsed beams of '^^Ti at the three bom¬ 
barding energies (300, 450 and 600 MeV) impinged on a 
metallic ^°Ca target (500pg/cm^ thickness) sandwiched 
between two very thin C foils (I5pg/cm^) to prevent 
prompt oxidation. Typical beam currents of 0.5-1 pnA 
were used. In order to study the GDR evolution with 
excitation energy and spin in molybdenum nuclei, a com¬ 
posite apparatus was used. A group of 8 BaF 2 scintilla¬ 
tors (Hector setup [1011211122]) for gamma rays, cov¬ 
ering backward laboratory angles, was coupled with a 
large acceptance detector for charged reaction products. 
Here we describe only this latter part of the setup, while 
details on the gamma array can be found in uniin]. 

The heavy products, in particular evaporation residues 
from the fusion reactions, were detected by an array of 
48 triple-phoswiches of the Fiasco setup [13| mounted 
in six matrices, each with eight detectors, in an axially 
symmetric configuration around the beam direction. The 
phoswiches featured two plastic layers (made of a fast 
180 pm, and a slower 5 mm scintillators) followed by a 
CsI(Tl) crystal 4 cm thick, thus presenting a wide dy¬ 
namic range for pulse shape analysis. With respect to 
their original use, these detectors were now equipped 
with digital electronics [ 23 ], purposely developed by the 
collaboration. The anode current pulses from the pho¬ 
tomultipliers were digitized and the relevant information 
(time mark, pulse shape and energy-related variables) 
was extracted. In practice, for each phoswich one ob¬ 
tains, besides a time mark for the time-of-flight (ToF), 
three energy-related variables (henceforth gateA, gateB 
and gateG); gateC corresponds to the light emitted by 


the GsI(Tl) with the longest time constant; gateB in¬ 
cludes the light of the second plastic layer and part of the 
light emitted by the GsI(Tl); gateA includes the fastest 
light components from all the three scintillator layers. 

The phoswiches covered the polar region from around 
5° to 25°; 32 of them were located in a wall configura¬ 
tion from 5° to 13° with a significant efficiency for ER 
detection while the remaining 16 scintillators were placed 
in two side arms, mainly to detect fission fragments. 
Thanks to the first fast plastic layer and to the large 
distance from the target ( 1 . 6 m), the phoswiches permit¬ 
ted velocity measurement of the ejectiles from Z = 1 
up to the ER with good time resolution (of the order 
of ns). The charge identification, only possible for the 
products punching-through the first thin layer, has been 
obtained via digital pulse shape methods from protons 
to Z close to the beam atomic number. The upper limit 
in charge depends on the particular telescope and on the 
beam energy, and it never exceeds Z = 18. Hydrogen iso¬ 
topes have been identified also by means of the fast-slow 
technique in GsI(Tl). Heavy fusion residues and heavy 
fission fragments, stopped in the first plastic layer, were 
detected without charge identification. It has to be noted 
that the time mark given by the phoswich telescopes by 
means of the digital technique (based on the numerical 
implementation of the Constant Fraction Discrimination 
method, used for the first time during this experiment) 
was progressively delayed when the kinetic energy of the 
particles increases and thus the energy deposited in the 
slow plastic layer and in CsI(Tl) grows. This systematic 
distortion (walk in the time mark detection) was associ¬ 
ated with a too small amount of digitized samples which 
were recorded for each phoswich signal and could not 
be corrected in the off line analysis. As a consequence, 
in the following, we discuss velocity spectra for heavy 
fragments which are stopped in the first layer, while for 
light products, entering also in the successive layers of 
the telescopes, we are able to rely only on the particle 
multiplicities; in fact we judge the velocity (and thus the 
energy) not trustworthy. 

The forward chamber of Garfield [2S], covering 
polar angles between 29.5° and 83° is able to detect 
light charged particles (LCPs) and intermediate mass 
fragments (IMFs). Garfield is based on AE{gas) — 
E(CsI(Tl)) modules, equipped with digital electronics 
and its features and performances are extensively de¬ 
scribed elsewhere [2S]. Very briefly, Garfield consists 
of 192 AE{gas) — i5(GsI(Tl)) telescopes. The gas volume 
is unique for all the modules and it is filled with flowing 
GF 4 at a pressure around 50 mbar. The collecting anodes 
of the AE stage, based on metal/glass microstrip tech¬ 
nology, provide a moderate gas multiplication. While 
preserving the linear response with deposited energy, the 
internal gain allows for detection and identification from 
light charged particles to intermediate mass fragments 
(typically Z « 10 - 12). The 192 AE{gas) - F;(GsI(T 1)) 
telescopes are organized in four polar rings correspond¬ 
ing to four GsI(Tl) shapes; each ring consists of 48 ef- 


4 


fective sectors, with an azimuthal granularity of 7.5°. In 
this experiment the mechanical constraints, related to the 
coupling of the BaF 2 with Garfield, imposed a change 
in the target holder. As a result the backward ring of 
Garfield (from 67° to 83°) was affected by the target 
shadowing and has been discarded in the present analy¬ 
sis. 

As already said, for each telescope the AE — E method 
allows to identify charged particles from Z = 2 to 
Z « 10—12, measuring also their energy. Moreover, 
light charged particles {Z = 1,2) were isotopically re¬ 
solved via pulse shape analysis in the CsI(Tl) crystals 
[2S] . The energy calibration of CsI(Tl) has been carried 
out by exploiting the many reference points collected in 
various measurements of elastic diffusion with low energy 
beams. The overall uncertainty on the energy measured 
by the CsI(Tl) is « 3-4%. For kinematic reasons, heavy 
residues cannot reach the angles of the Garfield setup; 
as a consequence, Garfield was able to identify (and 
measure the energy of) all the charged products hitting 
its active surface. 

Finally, a small plastic scintillator was located at « 
2 °, well below the laboratory grazing angle for all the 
investigated reactions (the most forward one is i?graz = 
6.9° for the 600MeV reaction). It covered a small solid 
angle (3.6 • 10“® sr) and has been used to count elastically 
scattered ions for absolute cross section normalization. 


III. EXPERIMENTAL RESULTS 

As stated in the introduction, the topic of this work 
consists in the investigation of the decay of the excited 
compound nucleus formed in the “^^Ti -|- "^°Ca reaction at 
three beam energies and the comparison with the predic¬ 
tion of the statistical code. The main reaction parame¬ 
ters for the investigated systems are reported in Table |l] 
In the columns labelled estimations 

of the maximum angular momentum for the existence 
of the compound nucleus are reported. The former is 
calculated in the framework of the rotating liquid drop 
model (RLDM) of [57], while the latter is obtained ac¬ 
cording to the modified rotating liquid drop model of 
Sierk |5S|. At 300 MeV, where the grazing angular mo¬ 
mentum Zgraz (obtained from [29]) does not exceed by a 
large amount ^^eTk/RLDM’ complete fusion represents 
the dominant reaction channel. The formed compound 
nucleus has a broad spin distribution m and it decays ei¬ 
ther by evaporation (fusion-evaporation (FE) process) or 
by fission (fusion-fission (FF) process, possibly followed 
by the evaporative decay from the Hssion fragments). 
The two fission products can be of similar mass (sym¬ 
metric fission) or may be formed by a lighter and a heav¬ 
ier fragment (asymmetric fission). The FF channel be¬ 
comes more and more significant with respect to the FE 
process when the angular momentum of the compound 
nucleus increases, up to the ^^eTk/RLDM value, beyond 
which the system becomes unstable against fission, be¬ 


cause the fission barrier vanishes. Beyond ^lierk/RLDM 
fusion is prevented and other reaction mechanisms take 
place; in particular the deep inelastic process is the most 
important, but also the quasi fusion / quasi fission mode 
(in which the system does not go through a completely 
equilibrated fused phase) can be present. Those mech¬ 
anisms, which are supposed to increase with beam en¬ 
ergy, produce two heavy residues in the exit channel (for 
the deep inelastic process, the quasi projectile QP and 
the quasi target QT), similar to symmetric fission. The 
contribution of such binary channels represents a back¬ 
ground for the true FF events which can be barely disen¬ 
tangled, because the populated phase space regions are 
partly overlapping. These mechanisms affect mainly the 
symmetric fission case because the two reaction partners 
are quite similar in terms of mass. 


A. Selection of fusion events 

In order to properly select FE and FF events, a first 
pre-sorting was done at the trigger level. Indeed, the 
main trigger was the coincidence of a hit in a phoswich 
and a hit in Garfield. Considering the phoswich and 
Gareield acceptances, this condition strongly selects 
FE events. In order to better suppress undesired events 
(i.e. elastic scattering and binary peripheral reactions) 
only detector signals above the noise level but below a 
given threshold (corresponding to gateA values around 
2000 ch, thick horizontal line in Eig. [^ were allowed to 
start the data acquisition. In Fig.j^ajwe show a typical 
gateA vs. ToE correlation, which is the main criterion for 
the off-line event selection, while in part (b) of the pic¬ 
ture the gateA vs. gateB correlation is presented. The 
plots refer to the reaction at 600 MeV. In the gateA vs. 
ToF correlation three different regions can be identified: 
ions stopped in the first thin plastic scintillator layer (hy¬ 
perbolic branch at ToF > 40 ns), for which the charge 
identification is not possible; IMFs and LCPs stopped in 
the second scintillator layer (spots at different gateA val¬ 
ues), for which the charge identification is possible and 
finally LCPs (only protons for the reactions at 300 MeV) 
entering the CsI(Tl) layer (sudden increase of the gateA 
value). Ions stopped in the first scintillator layer produce 
the almost vertical and unresolved branch in the gateA 
vs. gateB correlation, while those punching through the 
fast plastic collapse on the almost horizontal lines, corre¬ 
sponding to the various Z, visible in the picture. LPCs 
entering the CsI(Tl) give rise to the discontinuity in the 
lower part of the gateA vs. gateB correlation. 

FE events are selected by requiring that only one frag¬ 
ment in the phoswich wall be inside the inner dashed-line 
areal gate drawn on the gateA vs. ToF correlation in the 
region of unresolved heavy products. In the same fig¬ 
ure one can clearly see the yield drop above the trigger 
upper level. The island of events at the highest gateA 
values corresponds to residual elastic scattered ions and 
has been used for ToF calibration. 
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Table I. Main reaction parameters. Eh is the beam energy; Ecm is the available energy in the center of mass, after removing 
the energy lost in half target; i9graz and Zgraz are, respectively, the grazing angle and the grazing angular momentum, calculated 
according to |29| : critical angular momentum beyond which the fission barrier vanishes, calculated according to 

|27| : Zgf,7k i® critical angular momentum beyond which the fission barrier vanishes, calculated according to [28) : i?cN is 
the excitation energy of the compound nucleus, calculated according to the formula Eq^ = Ecm + Q, where Q is the Q-value 
for fusion; E^j^jA is the excitation energy per nucleon of the compound nucleus; ctr is the reaction cross section calculated 
according to |29 |. 
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Figure 1. (Color online) (a): correlation between gateA (see text) and the time of flight for a phoswich telescope at ~ 8.3°; 
(b): correlation between gateA and gateB for the same phoswich telescope. Areal gates for particle identification are shown in 
the left panel and are explained in the text. Data refer to the reaction at 600 MeV. 


The fission of nuclei with A « 100, below the Businaro- 
Gallone point, mainly produces two fragments with dif¬ 
ferent sizes. Therefore we expect that in our system the 
asymmetric splits prevail, as reported in |30j for molyb¬ 
denum isotopes. 

Candidates for symmetric FF events have been selected 
requiring the coincidence of two heavy products in the 
phoswich wall falling inside the outer full-line areal gate 
drawn on the left side of Fig. instead candidates for 
asymmetric FF events have been selected by the coinci¬ 
dence of a heavy product in the phoswich wall (inside 
the same full-line areal gate) and an IMF in whatever 
detector section. Moreover, candidate fission fragments 
must have their center of mass velocity collinear and 
their relative velocity compatible with the fission sys- 
tematics ED ED- For kinematical reasons, two heavy 
fragments (symmetric fission) cannot reach Garfield, 
while only the lighter fragment from asymmetric splits 
can be detected in Garfield. Whatever the emission 
angle, the IMF is the only partner of the fission process 
which can be identified in charge. Of course, asymmetric 


fission events in which the light partner is missed pro¬ 
duce a background of incomplete events for the selected 
FE event set. This background has been estimated, as 
described later. 


B. Use of the simulation 

We first consider the comparison of the experimental 
observables with Gemini-] — h in the reaction at 300 MeV, 
where the assumption of the formation of the GN fol¬ 
lowing complete fusion is most reliable; in particular, we 
compared experimental energy spectra and multiplicities 
of protons and a-particles with the prediction of the code 
for the decay chain of a ®®Mo source with the proper ex¬ 
citation energy (1.4 MeV/nucleon) and a triangular spin 
distribution up to the Zcrit value. In order to perform 
a correct comparison, the simulated events were filtered 
with a software replica of our setup. The main code pa¬ 
rameters were tuned in order to obtain the best overall 
agreement between experimental results and simulated 
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Figure 2. (Color online) (a): ER lab velocity distribution at the three bombarding energies (black: 300MeV, red: 450MeV, 
blue: 600MeV); symbols are the experimental data, while continuous lines correspond to GeminiH— h simulation with RLDM 
yrast, RLDM fission barrier, w = 1.0fm, Td = Ozs (see text). Part (b): ER lab angular distribution at 300MeV (black) and 
600MeV (blue); symbols are experimental data, while continuous lines correspond to GeminiH— h simulation. In both pictures 
dotted lines correspond to simulated fission events in which the light partner was lost and the heavy one was erroneously 
identified as FE. All spectra are normalized to their integral, except for dotted lines (for which the same scaling factor of 
continuous curves has been used). 


data. In particular, referring to the parameters intro¬ 
duced in [8], we tried to tune the level density param¬ 
eter deff (acting, in particular, on the parameter fcp of 
equation (15) of reference [5]), the spread of the LCP 
Coulomb barrier taken into account by a parameter w in 
the transmission coefficients, the yrast parametrization 
Ey, the time delay for fission and the parametriza¬ 
tion of the fission barrier Bf as a function of the total 
angular momentum J. According to [8], the level density 
parameter rules the slope of the exponential tail of the 
energy spectra while w controls their shape close to the 
Coulomb barrier. Ey is particularly significant for light 
nuclei with small moment of inertia; it affects the emis¬ 
sion of a particles which remove large amount of spin 
and has slight influence on protons. The ru parameter 
and the shape of Bf can influence the occurrence of the 
fission process and its competition with particle decay. 
The model parameters have been tuned mainly looking 
at LCP spectra for FE events. 

Once a reasonable parameter set has been fixed, we 
can use simulated events to reliably evaluate the detec¬ 
tor efficiency and thus estimate the 47r-corrected particle 
multiplicities and the absolute cross sections at all the 
measured beam energies. 


C. Fusion-Evaporation and Fusion-Fission 

According to the GeminiH— h simulation, the detection 
efficiency for the evaporation residues (ER) in FE events 
has been evaluated to be within 10 % and 13 % (increas¬ 
ing with the beam energy) with the parametrization giv¬ 
ing the best agreement between experimental and simu¬ 
lated LCP spectra (see in the following). In Fig. [^a) the 
experimental lab velocity distribution for the ER (full 


symbols) is compared with the simulated one (filtered 
with a software replica of our setup, continuous lines) for 
the three different beam energies; spectra are normalized 
to their integral. Simulated and experimental data are 
generally well matched, with the partial exception of the 
600 MeV case; as a consequence, we can argue that a 
residual pollution of other reaction mechanisms (not in¬ 
cluded in the simulation) is present for the highest beam 
energy. On part (b) of the same figure the experimental 
lab angular distribution of the ER (full symbols) for the 
reactions at 300 MeV and 600 MeV compared with the 
GeminiH— h simulation (continuous lines) is plotted. A 
very good agreement is obtained for this observable. In 
both sides of Fig. [^dotted lines correspond to a fission 
background of ER events where one of the fission prod¬ 
ucts (the lighter one) was lost and the other one fell inside 
the ER detection gate (dashed contour in Fig.[^. In this 
case the scaling factor is the same used for the continuous 
curves. These events show an almost flat angular distri¬ 
bution and they constitute a background of incompletely 
detected events, present both in simulated and experi¬ 
mental data. According to the GeminiH— h calculation, 
their amount is about 9 % of ER events at 300 MeV, it 
rises to about 34% of ER at 450 MeV and it is about 
50 % of the total ER at the highest beam energy. 

Protons and a-particles represent the dominant con¬ 
tribution in the decay chain of ER in FE events. The 
experimental center of mass energy spectra of protons 
and a-particles at all the investigated beam energies are 
presented in part (a) and (b) of Fig. respectively; the 
spectra are normalized to the number of ER and refer 
to particles detected in the Garfield ring covering the 
polar angle between 41° and 52°. 

As expected for the decay of a hot nuclear source, the 
spectra have Maxwellian shape, with the apparent tern- 
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Figure 3. (Color online) (a): experimental center of mass energy spectra for protons in the reaction at 300MeV (black 
squares), at 450MeV (red up triangles) and at 600MeV (blue down triangles), detected in the Garfield ring with polar angles 
ranging from 41° to 52°. The spectra are normalized to the number of evaporation residues, (b): same for a-particles. 


perature (given by the inverse slope of the high energy 
tail) increasing with the beam energy, as the excitation 
energy of the CN increases. 

A comparison of the experimental center of mass ki¬ 
netic energy spectra for light particles (black squares) 
with the prediction of the GeminiH — h code (red curve), 
run with standard values of the parameters (those re¬ 
ported in 0 ) and filtered with a software replica of the 
setup, is plotted in Fig. for the reaction at 300 MeV 
for particles detected in the Garfield ring covering the 
polar angle between 41° and 52°. Side (a) concerns pro¬ 
tons, while on side (b) a-particle spectra are presented. 
The spectra are normalized to their integrals, in order to 
put into evidence differences in their shapes. 

From these plots a good agreement emerges between 
the model and experimental data in the case of protons, 
while the a-particle spectra are strongly different: in fact 
the experimental Coulomb barrier is higher than the sim¬ 
ulated one. 

In order to improve the model agreement in the case 
of a-particles, we investigated the effect of some model 
parameters, chosen among those mainly influencing the 
shape of the kinetic energy spectra of evaporated parti¬ 
cles as reported in [5]. 

The effect of the level density parameter was investi¬ 
gated varying in the range 6-10 MeV. According to 
[5] the standard ko value for ®®Mo CN is 7.3 MeV. Please 
note that according to equation (15) of [8], since the n 
parameter is close to 0.1 for A = 88, for our system the ef¬ 
fective level density parameter a is almost independent of 
the excitation energy and approximately equal to d = ^ 
m- We found that for protons the standard value ko of 
[5] is fine, while no appreciable improvement of the agree¬ 
ment between simulated and experimental spectra could 
be obtained for a-particles in the investigated range. As 
a consequence, the standard value was kept. 

In Gemini-] — h the effect of the source deformation due 
to thermally induced shape fluctuations is taken into 
account by averaging the transmission coefficients over 


three different values calculated with three different ra¬ 
dius parameters of the nuclear potential (Rq, Rq — Sv and 
Rq -I- Sr, equation (9) of reference [S]), where Sr = wy/T 
and T is the nuclear temperature of the daughter nucleus. 
The spread of the transmission coefficient mainly influ¬ 
ences the shape of TCP energy spectra in the Coulomb 
barrier region and its effect is evident for a spectra, while 
it is almost negligible for protons, due to their lower 
Coulomb barrier. The standard Gemini-]— h prescription 
is w = 1.0fm; however, in [5] a case (^°®Cd, [13]) is pre¬ 
sented where ic = 0 fm (corresponding to a single barrier, 
i.e. a spherical nucleus) fits better. Note that the author 
of [8] doesn’t exclude that this discrepancy is due to a 
possible contamination from reaction mechanisms differ¬ 
ent from FE. In any case, on the basis of this evidence, we 
investigated the effect of w on the simulated spectra in 
the range w = 0.0-1.5 fm. Proton spectra are not affected 
by the variation of this parameter, while for a-particles 
a slight improvement (consisting in a shift of the barrier 
towards the experimental value) has been obtained by 
decreasing the w parameter. 

The other parameter able to influence the shape of en¬ 
ergy spectra mainly for a-particles and mainly for light 
systems {A < 150) is the shape of the yrast line; in [S] 
a shape that follows the Sierk parametrization |28] up to 
a threshold angular momentum and then increases lin¬ 
early is proposed for light systems (as our ®®Mo), obtain¬ 
ing very good agreement with many experimental spec¬ 
tra. However, it has to be noted that the explored spin 
range in [5] (up to 40-50 h) is below the angular momenta 
reachable by our system (see Table]!]) and thus this recipe 
may not work properly in our case. Another shape of the 
yrast line is delivered by the RLDM of m- 

The two different yrast recipes ([5] and RLDM [57]) for 
our system are shown in Fig. ]^ as solid curves. The two 
curves are very similar up to J ~ 35 h; above this value 
the RLDM yrast (red thick curve) tends to be similar 
to that of a spherical nucleus (not shown). In the same 
Fig. ]^we compare also two recipes for the fission barrier. 
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Figure 4. (Color online) (a): center of mass energy spectra of protons in the reaction at 300MeV for experimental data 
(black squares) and Gemini++ simulation (red curve) with standard parameters (from [8]). The spectra are normalized to 
their integrals. Data refer to particles detected in the Garfield ring covering the polar angles between 41° and 52°. (b): same 
for a-particles. 


Yrast and fission barriers 



Figure 5. (Color online) Solid curves: yrast line parametriza- 
tions for ®®Mo nuclei as a function of their angular momen¬ 
tum. Black thin curve: linearized Sierk yrast (prescription of 
m)- Red thick curve: Rotating Liquid Drop Model yrast Ha- 
Dotted curves: Fission barrier parametrizations for ®®Mo nu¬ 
clei as a function of their angular momentum. Black thin 
curve: Sierk fission barrier [^. Red thick curve: RLDM 
fission barrier Ba¬ 
the Sierk prescription [55] (black thin dotted curve) and 
the RLDM (red thick dotted curve). At a fixed an¬ 
gular momentum the RLDM barrier is higher than the 
Sierk one, thus reducing the fission probability for the 
ER. This has influence not only on the shape of the par¬ 
ticle energy spectra but also on their multiplicities. 

Proton spectra are almost independent of the adopted 
parametrization, while for a-particles the best results 
are obtained when the RLDM barrier is coupled to the 
RLDM yrast, as shown in Fig. (blue thick dash-dotted 
curve). We mention that this particular Gemini-] — j- pa¬ 
rameter set corresponds to the choice made in HU. 

The agreement for the energy spectra of a-particles can 
be further improved if we couple the RLDM barrier and 



Figure 6. (Color online) Center of mass energy spectra 
of Q-particles in the reaction at 300 MeV for experimental 
data (black squares) and GeminiH— h simulations with differ¬ 
ent prescriptions for the fission barrier and the yrast curve. 
The standard value (Sierk fission barrier and Linearized Sierk 
yrast) is the red continuous thick curve; the green dotted 
curve corresponds to standard yrast (Linearized Sierk yrast) 
and RLDM fission barrier; the thin blue curve corresponds 
to standard fission barrier (Sierk fission barrier) and RLDM 
yrast; the blue dash-dotted curve corresponds to RLDM fis¬ 
sion barrier and RLDM yrast. Spectra are normalized to their 
integrals. Data refer to particles detected in the Garfield 
ring covering the polar angles between 41° and 52°. 


RLDM yrast prescription with the choice of w = Ofm, 
i.e. if we switch off thermal barrier fluctuations, moving 
towards a spherical nucleus, as it is shown in Fig. blue 
curve; in particular, the better agreement obtained in the 
barrier zone (zoomed in the inset) has to be noted. 

The time delay for fission ru has no influence on the 
particle energy spectra (Fig. comparison between red 
continuous and red dotted curve), while it affects the fis¬ 
sion probability and the angular distribution of the emit- 
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Figure 7. (Color online) Center of mass energy spectra of 
a-particles in the reaction at 300 MeV for experimental data 
(black squares) and Gemini++ simulations with different pre¬ 
scriptions. Blue curve: RLDM fission barrier, RLDM yrast 
line, w = 1.0 fm and Td = Ozs. Red curve: RLDM fission 
barrier, RLDM yrast line, w = Ofm and Td = Ozs. Red dot¬ 
ted curve: RLDM fission barrier, RLDM yrast line, w = Ofm 
and Td = 10 zs. In the inset the barrier region is zoomed with 
a linear ordinate scale. Spectra are normalized to their inte¬ 
grals. Data refer to particles detected in the Garfield ring 
covering the polar angles between 41° and 52°. 


ted particles (as it will be discussed in Fig. 11). 

The experimental center of mass energy spectra of pro¬ 
tons (panels (a), (b), (c)) and a-particles (panels (d), (e), 
(f)) for three angular rings of Garfield are compared in 
Fig.j^with the corresponding simulated distributions ob¬ 
tained with the “best” set of parameters so far discussed 
{ko = 7.3 MeV, RLDM fission barrier, RLDM yrast line, 
w = Ofm, Td = 10zs). The agreement between exper¬ 
imental data and simulation is reasonable at all angles 
accessible with Garfield. 


Calculations have been extended to the two higher 
beam energies assuming the parameter set tuned for the 
300 MeV reaction. 

The center of mass energy spectra of protons (panels 
(a), (b), (c)) and a-particles (panels (d), (e), (f)) for the 
reactions at 450 MeV and 600 MeV are displayed in Fig.[^ 
and Fig. respectively. Also in this case a reasonable 
agreement at all the explored lab angles is observed be¬ 
tween simulated and experimental data. 

Laboratory angular distributions for protons (panels 
(a), (b), (c)) and a-particles (panels (d), (e), (f)), normal¬ 
ized to the number of ER, are shown in Fig.[^for all the 
investigated reactions (panels (a) and (d): 300 MeV; (b) 
and (e): 450 MeV; (c) and (f): 600 MeV). Experimental 
data are plotted as black full circles, while different Gem¬ 
ini-]— j- calculations, all assuming RLDM yrast line and 
RLDM fission barrier, correspond to different colors. The 
horizontal bar in the Garfield region ('d > 30°) indi¬ 
cates the angular coverage of each ring. The vertical error 
bar for experimental data represent estimated systematic 
errors due to the evaluation of the efficiency in the in¬ 


dividual detection cells; statistical errors are negligible. 
The main observations are as follows. Concerning proton 
emission, the data are overall well described by the statis¬ 
tical model at the three energies, whatever the input pa¬ 
rameter choice. There are minor discrepancies which are 
different for the three bombarding energies; at 300 MeV 
the model underestimates the proton yields beyond 30°, 
while the opposite happens at 600 MeV. However the 
measured proton emission appears to be essentially com¬ 
patible with the evaporation from ®®Mo compound nu¬ 
clei. The most striking difference between experiment 
and model is for a emission, as the measured distribution 
becomes more forward peaked than the model predictions 
when the energy increases, almost independently of the 
fine tuning of the parameters. There are of course some 
changes if some parameters are modified. For example, 
the best set from the point of view of the shape of the 
energy spectra (w = Ofm, Td = lOzs, red squares) gives 
the worst agreement for the forward emission. The choice 
w = 1.0fm (green up triangles and blue down triangles), 
corresponding to a slightly poorer reproduction of the a 
energy spectra (see Fig. [^, increases the forward focus¬ 
ing of a-particles, with a minor effect of the delay time 
for fission (rd = Ozs for blue down triangles, Td = 10zs 
for green up triangles). In any case, the experimental 
focusing of a-particles, increasing with the beam energy, 
cannot be properly reproduced by the model code. 

Concerning the origin of the observed yield discrep¬ 
ancy for a-particles we have tested several hypotheses. 
For example, the possibility that the "^^Ca target is oxi¬ 
dized (and, as a consequence, causes a spurious emission 
of a particles) has been excluded by means of a chemical 
analysis on the sample. A possible source of a contami¬ 
nation comes from the C layers embedding the ^°Ca; the 
detected events can contain a contribution from the com¬ 
plete fusion of “^^Ti and This hypothesis has been 
rejected too, because of the very different center of mass 
velocity of the ^®Ti -j- system with respect to ®®Mo 
(they differ by about 12 mm/ns at 600 MeV and about 
9mm/ns at 300MeV): events coming from the complete 
fusion of 4®Ti and are well outside the selected ER 
gate (Fig. 0- The fission background of ER events (not 
negligible, mainly at 600 MeV) is already included in the 
Gemini-]— ]- simulation; any attempt to reduce this spu¬ 
rious contribution by means of a strict cut on the most 
forward emitted ER has little effect since Gemini-]— h 
shows that the angular distribution of the background 
fission fragment is rather broad and flat. In any case, 
as shown by Fig. for a-particles at 600 MeV, we have 
verified that, according to Gemini-]--)-, particles emitted 
by a fission fragment erroneously identihed as ER are 
only weakly more forward focused (full black dots) with 
respect to the case of true ER (red squares); as a con¬ 
sequence, they cannot explain the observed discrepancy 
between simulation and experimental data. 

Concerning Deep Inelastic Collisions, not included in 
the adopted model, we note that for kinematical reasons 
only a small tail of QP/QT from very dissipative colli- 
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Figure 8. (Color online) Center of mass energy spectra for LCPs detected at different lab angles for the reaction at 300 MeV. 
Experimental data correspond to black squares. Simulated data have been obtained with ko = 7.3 MeV, RLDM fission barrier, 
RLDM yrast line, w = Ofm, Td = 10zs. Spectra are normalized to their integral, (a), (d): 53.0° < diab < 66.0°; (b), (e): 
41.0° < i9iab < 52.0°; (c), (f): 29.5° < diab < 40.0°. Panels (a), (b), (c) correspond to protons, while panels (d), (e), (f) 
correspond to a-particles. 




450MeV p 47° 




450MeV p 35° 




Figure 9. (Color online) The same as Fig. [^but for the reaction at 450MeV. 


sions could fall inside in the ER identification gate. As 
a consequence the possible contamination should be neg¬ 
ligible, although a simulation of the DIG process based 
on a phenomenological parametrization of this mecha¬ 
nism indicates a preferential emission of a-particles in the 
forward direction at 600 MeV. Other processes such as 
QuasiFusion/QuasiFission are not included in the model 
and could contaminate in a minor way the experimental 
sample of ER events. 

With the above limitations, we can use GeminiH— h 
data in order to estimate the efficiency for LCP detec¬ 
tion and to correct the measured yields. Various slightly 
different correction factors have been deduced according 


to the different adopted parametrizations of Gemini-]—]-• 
The resulting LCPs 47r-multiplicities obtained by aver¬ 
aging over the values obtained with different parameter 
sets are reported in Fig. for protons, a-particles and 
deuterons, exploiting the isotopic resolution for hydrogen 
isotopes both for Gareield and the phoswich telescopes; 
tritons are too weak to be reliably evaluated. 

Vertical error bars correspond to the standard devia¬ 
tion of the different results. Each point is also entitled to 
a statistical error of the order of 8 % for protons (circles), 
20% for deuterons (triangles) and 40% for a-particles 
(squares). The plot shows that particle multiplicities 
slightly increase with the beam energy, with a weaker 
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Figure 10. (Color online) The same as Fig. but for the reaction at 600 MeV. 
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Figure 11. (Color online) Differential cross section as a function of the lab polar angle. Panels (a), (b), (c): protons. Panels 
(d), (e), (f): a-particles. (a) and (d): 300 MeV; (b) and (e): 450 MeV; (c) and (f): 600 MeV. Full black circles: experimental 
data. Red squares: Gemini++ with RLDM yrast line, RLDM fission barrier, w = Ofm, Td = 10 zs (best choice for energy 
spectra). Blue down triangles: GeminiH— h with RLDM yrast line, RLDM fission barrier, w = 1.0fm, Td = Ozs. Green up 
triangles: GeminiH— f with RLDM yrast line, RLDM fission barrier, w = 1.0 fm, Td = 10 zs. 


trend for protons. The observed forward focusing of a- 
particles increasing with bombarding energy suggests the 
onset of pre-equilibrium processes. Therefore, we tried 
to estimate a multiplicities separately from Garfield 
data only (extrapolating them to the full solid angle), 
where angular distributions and yields are quite well re¬ 
produced by the statistical model. The same choice can 
be applied also to protons, although no clear evidence of 
a model failure with increasing energy exists (Fig. 11 top 
part). The obtained results are shown as open symbols 
in the same Fig. For a-particles the difference be¬ 
tween the two data sets represents an upper limit for pre¬ 
equilibrium emission that, as expected, increases with the 


beam energy, or, in any case, for emission sources not in¬ 
cluded in the simulation. On the contrary, for protons 
there is no significant discrepancy at all energies; even 
at 600 MeV the two extracted proton multiplicities are 
compatible within the errors. For a-particles the possi¬ 
ble pre-equilibrium contribution starts at 450 MeV and 
corresponds to about 0.5 particles per event (with re¬ 
spect to a total number of 2.8 a-particles per event) in 
the worst case (600 MeV), thus justifying the assumption 
of negligible pre-equilibrium effects done in m for the 
evaluation of the GDR strength. 

By means of the fast plastic scintillator located below 
the grazing angle for all reactions (see section 0 it has 
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Figure 12. (Color online) Simulated angular distribution for 
a-particles at 600 MeV; Gemini++ data with RLDM yrast 
line, RLDM fission barrier, w = 1.0 fm and ra = Ozs. Red 
squares correspond to events in which a true ER is selected, 
while black dots correspond to incomplete events in which a 
fission fragment produces a background of the fusion channel. 
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Figure 13. Estimated experimental particle multiplicities in 
Itt as a function of the beam energy in FE events. Circles; 
protons; triangles: deuterons; squares: a-particles. Full sym¬ 
bols refer to data from Garfield and phoswich detectors. 
Open symbols give the particle multiplicities extracted from 
Garfield data only. 


been possible to estimate absolute cross sections. For FE 
events, the obtained values are reported in Table [III in the 
row labelled ctfe; they are the average among the various 
estimated ctfe values obtained running Gemini-]— j- with 
different sets of parameters (all of them with RLDM yrast 
line and RLDM fission barrier). The error is the standard 
deviation of these values. It is interesting to note that, as 
shown in the last row of the table, the FE cross section is 
almost one half of the reaction cross section (last column 
of Table |T]) at 300 MeV, while it decreases to about 19 % 
at 600 MeV. 

If we plot the normalized FE cross section ctfe/cr as 
a function of the available c.m. energy per nucleon as 
proposed in the recent systematics for fusion evaporation 


Table II. Experimental absolute cross sections, ctfe is the 
fusion-evaporation cross section; ctff is the fusion-fission cross 
section, ctf/ctr is the ratio between the total fusion cross sec¬ 
tion (FF plus FE) and the total reaction cross section, taken 
from the last column of Table |I] ctfe/ctr is the ratio between 
the FE cross section and the total reaction cross section. 



300 MeV 

450 MeV 

600 MeV 

CTpE 

(893 ± 109) mb 

(545 ± 45) mb 

(459 ± 115) mb 

CTFF 

(115 ±3) mb 

(266 ± 37) mb 

(417 ± 114) mb 

CTf/ctr 

0.54 ±0.06 

0.358 ±0.004 

0.36 ±0.07 

CTfe/ctr 

0.48 ±0.06 

0.24 ±0.02 

0.19 ±0.05 


[9], we can see (Fig. 14) that our data (red dots) fairly 
agree with the reported prescription. 



Figure 14. (Golor online) FE cross section normalized to 
the reaction cross section as a function of the Ecm/nucleon. 
Red dots correspond to the experimental values found in this 
work, compared with the systematics of [9]. The figure is 
adapted from Fig. 2 in [9]. 


Finally, we briefly discuss the fission case in order to 
extract the absolute cross section for this channel. The 
fission process is mainly associated with the highest an¬ 
gular momenta, as it is evident from Fig. dotted lines. 
In the adopted parametrization (RLDM fission barrier) 
the barrier is around 50 MeV at low J (as calculated also 
in [33] ) and decreases smoothly, until vanishes around 
79 h. In the experimental data the symmetric fission is 
weak, the measured fraction varying from 10 % to 30 % of 
the total fission events (without geometrical correction), 
depending on the bombarding energy. The absolute fis¬ 
sion cross section ctff is reported in Table [III second row: 
its weight in the total fusion cross section smoothly in¬ 
creases from 300 MeV to 600 MeV. At 600 MeV the FF 
cross section is of the same order of magnitude as the FE 
cross section. If we compare the total fusion cross sec¬ 
tion with the total reaction cross section ctf/ctr (third 
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line of Table 0 we can see that also at 300 MeV there 
is room for other reaction mechanisms (such as DIG): in 
fact fusion represents only 54 % of the total reaction cross 
section. 


IV. SUMMARY AND CONCLUSIONS 

In this work we presented some exclusive data con¬ 
cerning the system ^®Ti -|- "‘^Ca at three bombarding en¬ 
ergies (300, 450 and 600 MeV). Data have been collected 
by means of the Garfield setup coupled to a wall of 
phoswich telescopes of the Fiasco device. 

We have focused our analysis on the fusion channel 
whose collective dipolar excitation has been the subject 
of our recent investigation m- In particular we have 
shown that: 

• this channel represents more than 50 % of the total 
reaction cross section at 300 MeV and it goes down 
to about 40 % at 600 MeV. 

As expected, the formed CN decays through the evapo¬ 
ration of LCPs (fusion-evaporation event) or can undergo 
fission (fusion-fission event). After selecting FE events, 
we have constrained the main parameters of the statisti¬ 
cal code Gemini-|--|- comparing measured and simulated 
LCP distributions at 300 MeV where contributions from 
sources other than the CN should be negligible. We have 
seen that: 

• proton energy spectra are reproduced by Gem¬ 
ini-)— |- with the standard parameters; 

• in the case of a-particles it is necessary to use 
the RLDM parametrization both for the yrast line 
and the fission barrier, and to reduce the LCP 
Coulomb barrier fluctuations. This can indicate 
rather spherical excited nuclei as the sources of the 
detected LCP; 

• at the two highest beam energies there is an excess 
of Qf-particles at forward angles, which cannot be 
reproduced by the statistical model; 


• this excess, corresponding on the whole to about 
20 % of the total emitted a-particles at the 
highest energy (600MeV), may be due both 
to pre-equilibrium emission and to other pro¬ 
cesses (such as Deep Inelastic Collision, QuasiFis- 
sion/QuasiFusion), whose residual contamination 
cannot be definitely excluded due to the difficulty 
of a very clean fusion channel selection; 

• this excess of a emission can be taken as an upper 
limit for the pre-equilibrium emission of this sys¬ 
tem, thus justifying the assumption done in m of 
negligible pre-equilibrium. 

Finally we have extracted the absolute cross section 
for the FE and EE cases: 

• the FE cross sections are in reasonable agreement 
with the existing systematics reported in [5]; 

• the total fusion cross section is well below the total 
reaction cross section (estimated according to P5]L 
suggesting the presence of other processes, such as 
Die, sizably contributing at the highest angular 
momenta. 
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